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We analyze the impact of trap states in the oxide layer of a superconducting tunnel junctions, on 
the fluctuation of the Josephson critical-current, thus on coherence in superconducting qubits. Two 
mechanisms are usually considered: the current blockage due to repulsion at the occupied trap states, 
and the noise from electrons hopping across a trap. We extend previous studies of noninteracting 
traps to the case where the traps have onsite electron repulsion inside one ballistic channel. The 
repulsion not only allows the appropriate temperature dependence of l//-noise, but also the control 
of the coupling between the computational qubit and the spurious two-level systems inside the oxide 
dielectric. We use second-order perturbation theory, which allows to obtain analytical formulas for 
the interacting bound states and spectral weights, limited to small and intermediate repulsions. 
Remarkably, it still reproduces the main features of the model as identified from the numerical 
renormalization group. We present analytical formulations for the subgap bound-state energies, the 
singlet-doublet phase boundary, and the spectral weights. We show that interactions can reverse 
the supercurrent across the trap. We finally work out the spectrum of junction resonators for qubits 
in the presence of onsite repulsive electrons and analyze its dependence on microscopic parameters 
that may be controlled by fabrication. 



PACS numbers: 74.50.+r, 72.10.Fk 
I. INTRODUCTION 

Josephson junction circuits are promising candidates 
for the physical implementation of quantum comput- 
ing. They promise scalability and integrability from the 
very start. The grand challenge of realizing quan- 
tum computing in such a mesoscopic solid-state device is 
decoherence^. In the first decade of experimentation in 
this area, decoherence from the electromagnetic environ- 
ment has been understood and exquisitely controllecP^. 
Over the last years, it has become increasingly clear that 
the limits in coherence times are set by properties of the 
material comprising the Josephson circuit, in particular 
its surfaces and interlayer dielectrics PCU Q ne p r0 minent 
source of decoherence is critical-current noise. This noise 
has already been discussed in the context of optimizing 
high-resolution measurements in superconducting quan- 
tum interference devices (SQUIDs). Several models have 
been put forw ard, including local moments from glass 
physics^- 7 * 18 * 27 !, Kondo and other impurities^H22l and non- 
interacting Andreev trap states^J 3 ^. 

A typical Josephson junction, a thin, in most cases 
non-crystalline, insulator of oxide layer between two 
pure superconductor reservoirs, may have a fluctuating 
state that causes subsequent fluctuations in the critical- 
current. One can argue these fluctuations could be due to 
defects in the amorphous material of oxide layer. A single 
defect, i.e., a two-level system, can generate a Poissonian 
distribution of current noise, the so-called "random tele- 
graph noise", and a bath of these fluctuators with natu- 
ral energy level distribution may superpose into a non- 



Gaussian 1// current noise. The 1// noise has been ex - 
perimentally observed in the Josephson junction^ 13 * 33 * 34 !. 

There are two experimental results that illustrate the 
characteristics of two-level systems generating this noise. 
First, the 1/ f noise in a junction below the frequency 
10Hz displays a Independence on temperature.^ Sec- 
ond, a more precise measurement performed by Sim- 
monds et.al. 8 , and subsequently, many others, shows that 
a current-biased phase qubit during driving transitions 
between states |0) and |1) is partially blocked at some 
specific frequencies. This indicates the existence of mi- 
crowave resonators in the current noise spectrum on top 
of the 1// spectrum. These resonances were reported 
in Ref. to be highly suppressed by 80% if a single- 
crystal AI2O3 replaces the amorphous tunnel barriers. 
Nevertheless the nature of this noise and the reason for 
its enhancement in crystalline dielectrics is left unknown. 

Faoro and Ioffe in Ref. [30] hypothesized that the In- 
dependence of the 1/ / -noise sector is satisfied if 1) the 
two-level systems are shallow Kondo traps, and 2) the 
majority of these traps have bound states very close to 
the Fermi surface. Inside a tunneling junction there are 
typically a few O2 impurities left from the oxidization. 
These Kondo traps are impurities with a singly occupied 
level that carry a free spin. These magnetic impurities 
can generate localized bound states below the supercon- 
ducting gap. Two potential candidates for generating 
noise at very low energy (/ < 10 Hz) are (1) slow elec- 
trons and (2) nuclear spins; however due to the universal 
1// noise persistence only up to 10Hz, which is two or- 
ders of magnitude smaller than the typical energy scale 
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associated with nuclear spins energiesJ^SHsZl the nuclear 
spins cannot play a cruc ial rule in generating this noise. 
Recent experiment s ' 38 ! 39 ! as well as a density functional 
calculation 40 confirmed the nature of localized electronic 
levels to act as electronic spin-oriented traps. 

Constantin and Yu in Ref. |4T] studied the critical- 
current noise from the interaction of conduction electrons 
with a magnetic dipole localized inside the Josephson 
junction and showed that it is suppressed as the mag- 
netic dipole is located deeper inside the junction. A 
deeper study of this phenomenon through quantum field 
theory unveil more details for the cases of shallowly or 
deeply located magnetic impurities inside the junction. 
If the Kondo traps are coupled to only one of the su- 
perconducting reservoirs (i.e., located close to one of the 
superconductor-insulator interfaces) , they may affect the 
critical-current by (1) blocking out a part of supercur- 
rent and (2) generating microwave resonances from the 
coupling between conduction electron and the electron 
captured by the trapping center. In a single trap and 
for the case of non-interacting electrons this has been 
studied by de Sousaei alP^ where microresonators ap- 
pear at frequencies twice as the Andreev bound-state en- 
ergies. However, that study suffers from the absence of 
Kondo temperature as the electrons do not interact. In 
the first part of this paper, we extend the noise study to 
nonzero Kondo temperature. On the other hand, if the 
Kondo traps are coupled to both superconducting reser- 
voirs (i.e., located deep in the cener of junction) they af- 
fect the critical-current by (1) the quasiparticle poisoning 
of the Josephson current and (2) the current reversal flow. 
Faoro, Kitaev and Ioffe in Ref. |35]studied the quasiparti- 
cle poisoning contribution to charge and critical-current 
noise via a simple toy model. In the second part of this 
model, we study the current reversal flow effect. 

The inclusion of Coulomb interaction alters the An- 
dreev bound states to form bound states elsewhere inside 
the superconducting gap. As Faoro and Ioffe indicated,^ 
onsite interaction has a major impact on the 1// noise 
sector such that without it, the temperature dependence 
of the current noise cannot be achieved with a reason- 
able density of states. Given the possibility of low-lying 
energy levels near the Fermi sea in the presence of in- 
teraction, and assuming that the majority of traps are 
at low-lying energy levels, the interaction impacts the 
charge blockade phenomenon too. 

In this paper, we aim to tackle this problem by study- 
ing the critical-current noise generated by repulsive elec- 
trons. To achieve this goal, one needs precise solutions 
of the Anderson impurity model in the Josephson junc- 
tion. As will be described in more detail, the mean-field 
approximation to these problems often leads to unphys- 
ical result, thus one needs to resort to non-perturbative 
methods such as the numerical renormalization group 
(NRG)f^, the functional renormalization group (fRGjS. 
and density matrix renormalization group (DMRG)I 45 * 46 I 
. These methods have achieved the precise description of 
the competition between superconducting order and the 



effect of the impurity. Nevertheless, NRG and DMRG 
need heavy numerical analysis and fine-tuning of numer- 
ical parameters, and fRG suffers from many truncations 
necessary to numerically generate the appropriate self- 
energies. Unfortunately, there is no secure analytical so- 
lution that is as reliable as NRG and that accurately 
describes experiments. We study this problem analyti- 
cally within second-order perturbation theory. We de- 
rive the analytical weights associated with the subgap 
resonances in the interacting theory and interestingly re- 
produce some of the features of interacting theory known 
previously by numerical analysis of NRG and fRG. After 
a discussion of the range of validity of our approach, we 
highlight the noise spectrum due to blockade of current 
by occupied traps. 



A. 7r-junctions and interaction effects 

A Josephson junction is characterized by a supercur- 
rent I((f>) that is an antisymmetric, 27r-periodic function 
of the phase difference between the superconductors^. In 
the regime of weak and time-reversal-invariant coupling 
between superconducting reservoirs across the junction, 
e.g., in the nonmagnetic tunnel junctions, we have the dc 
Josephson current I = I c sin <j) with the phase difference 
of order parameter of the two superconducting reservoirs 
connecting via a junction and I c the critical-currenlPSl. 

The electron transport through a junction between two 
superconducting reservoirs is hea vily affected in the pres- 
ence of magnetic impurities EMS: This coupling breaks the 
time-reversal symmetry and leads to the appearance of 
states below the superconducting energy gap inside the 
junction. When the ground state is occupied by one elec- 
tron, another electron is repelled from it. This repulsion 
affects the current flow such that the conduction electrons 
may either screen the trapped electron or flip its spin. 
This is determined through the competition between elec- 
tronic interaction (characterized by Kondo temperature 
Tk) and superconductivity (characterized by the super- 
conducting gap A). The possibility of spin flip scattering, 
no matter if it is a dominant process or not, results in 
the modification of the bound-state energies such that it 
becomes possible that at a specific interaction strength 
the bound states cross the Fermi surface. 

In a metallic system the Kondo traps occupied by one 
electron spin at high temperature causes conduction elec- 
trons to scatter off inelastically in a spin flip process, 
while at low temperature, they are screened forming a 
bound singlet state leaving only elastic scattering. The 
Kondo temperature Tk is the temperature crossover be- 
tween the two different scatterings. However, in a super- 
conducting system more complication appears due to the 
competition between the Kondo temperature and the su- 
perconducting gap. When the repulsion between onsite 
electron and the conduction electrons is weak (i.e., large 
Kondo temperature, Tk 3> A) in the superconducting 
phase the magnetic moment is screened by the conduc- 
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tion electrons such that transport through the junction 
occurs without spin flip and the ground state of the sys- 
tem (superconductor + impurity) is singlet. In this case 
the electron at the magnetic impurity along with the con- 
duction electrons form a bound-state below the supercon- 
ducting gap at the impurity site. However, strong inter- 
actions (i.e.. low Kondo temperature. <C A) causes 
the conduction electrons to spin flip at the impurity site, 
where the ground state of the system is a doublet. This 
affects the dc current to flow opposite to the phase drop, 
i.e., dl((f>)/d(j)\ r p =0 < 0. This sometimes is referred to 
as a 7r-junction, as in this case the ground state at zero 
phase bias is characterized by 4> — tt. In general, the 
resulting current-phase relation is not a simple sinusoid. 
There is a "certain" Kondo temperature that serves as the 
border between spin flip phenomenon and the screening. 
Recently, using NRG method, this divider Kondo tem- 
perature was reported to be « 0.3A.^ 
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Figure 1. (Color online) (a) Energy view on a Josephson 
junction sandwiched between two superconducting reservoirs. 
The two thick (red color) bands are the subgap bound states 
caused by the presence of magnetic impurity in the junction, 
(b) Cross-section of the tunneling layer. Traps near the oxide 
surface (left) predominantly contribute to blocking noise — an 
occupied trap is charged and blocks off parts of the current. 
Traps in the center carry a contribution to the supercurrent 
that is influenced by interactions and can potentially lead to 
a 7r-junction. 

There have been various proposals and observations 
of 7r-junctions including those at high-temperature su- 
perconductors, non-equilibrium mesoscopic junctions, 
and magnetic junctions. Here, we are discussing an 
interaction-induced mechanism. When two supercon- 



ducting reservoirs are connected via a Josephson junc- 
tion whose width is sufficiently smaller than the coher- 
ence length, spin flip tunneling processes can cause su- 
percurrent reversaP2HS21 More recently van Dam and 
Nazarov observed this phenomena in a superconducting 
quantum dotF*H It was shown that at high Kondo tem- 
perature (i.e A <C Tx) the Josephson current through a 
Kondo impurity is non-sinusoidal. Ishizaka calculated the 
Josephson current numerically within the non-crossing 
approximation; it was reported that as a function of the 
Coulomb repulsion the supercurrent is suppresse d unti l 
at some critical value a sudden sign change occuriPSES , 
Within the numerical renormalization technique, it was 
reported that a 7r-junction is the preferred state in junc- 
tion when the impurity is singly occupied and the onsite 
Coulomb interaction is large^SHSD A generalization of 
the non-crossing approximation found a modification in 
this to 7r junction transition at the presence of bound 
states ^. Rec ently some other approaches including NRG 
methocPM] an d fRGP were taken to study this effect 
on a quantum dot and they predicted a 0-7r transition at 
large Coulomb repulsion. 

Figure Qa) illustrates bound-state energies (— e and 
U — e) inside the Josephson junction caused by magnetic 
impurity. In the particle hole symmetrically case the two 
states are located symmetric above and below the Fermi 
surface, i.e., e :~ U/2) The renormalized width of these 
levels is determined by the Kondo temperature. As indi- 
cated by Yoshida and OhashPSland confirmed by Bauer 
et. all^the Kondo temperature in the symmetric case is 
defined by T K = 0.182C/ v /8r/7rt7exp(-7r[//8r), with U 
being the Coulomb repulsion and T being the transition 
rate between superconductor and impurity. 

Figure Q shows transport inside a junction in the 
presence of a magnetic impurity. Let us first consider 
the first row of panels (al)-(a3). In weak-interaction 
limit, the impurity electron is screened by conduction 
electrons, therefore the other electrons are scattered only 
from the potential of this singlet and no spin exchange 
occurs at the impurity site. The panel al shows the bro- 
ken Cooper pair in the left reservoir due to proximity 
effect inside the insulator junction and either one of the 
two electrons can tunnel to the other reservoir [see panel 
(a3)] via the virtual state of panel (a2) without spin ex- 
change. The two electrons after consecutive tunneling 
reproduce the original Cooper pair now in the right side 
with parity conservation. In the strong repulsion regime 
(i.e., U > U c ), the impurity electron is unscreened and 
the conduction electrons caused by breaking Cooper pairs 
inside the junction exchange their spin with the impurity 
electron. As it is shown in the bottom row of Fig. ^ 
panels (bl)-(b3), due to the Pauli principle between the 
two outcome electrons the one with opposite spin first oc- 
cupies the impurity level and the impurity electron may 
tunnel out. Repeating the same process causes a spin flip 
on both electrons reproducing a Cooper pair in the other 
reservoir. In other words, the initial Cooper pair state 
(I tl) — I 4-t))/v2 after passing through this junction be- 
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comes e l7r (| t|) - | 4t))/V2- This simply results into a 
7r shift in the Josephson relation and a negative sign for 
the supercurrent. 



B. Occupation-blockade noise spectrum 

A more conventional junction noise mechanism is 
based on the occupation number of the trap state. The 
idea is that whenever the trap is occupied by one or 
two electrons. Coulomb repulsion originating from the 
localized trap will block off part of the junction area, 
hence influencing other channels. This requires the trap 
to be coupled to only one of the electrodes, whereas for 
the 7r-j unction scenario we require that the coupling of 
the trap to both electrodes is strong enough. Occuption 
noise has been stud ied early in the 1980sS3E3 anc j more 
recently^HanEUMMI. These mechanisms are illustrated 
in Fig. (|]). 

The charge-blockade decoherence determines how sen- 
sitive the microwave resonances, present on top of 1/f 
noise spectruirP^, are to the Coulomb repulsion. As a re- 
sult, both the energy of these resonances and their inten- 
sities are heavily affected in strong Coulomb repulsion. 

Our goal is to use an analytical approach to the prob- 
lem. It is well-known that the Kondo effect in normal 
metals cannot be described by a mean-field study be cause 
it is uncontrolled and underestimates correlations^^, 
however when this method was utilized by Yoshioka and 
OhashPSl on superconductors interestingly it was shown 
at a critical Coulomb interaction the current starts to 
flow backward, which means some — 7r transition is 
captured in mean field limits. Nevertheless, later on it 
was understood this is the artifact of an unphysical ef- 
fect because in the mean field approximation there is no 
well defined Kondo temperature that allows to control 
the phase transition from a non-magnetic phase into a 
magnetic one and this current reversal flow appears only 
because of the unphysical local moments^. 

We extend this study into the second-order pertur- 
bation in Coulomb repulsion and will show that the 
results are remarkably consistent wit h the recent non- 
perturbative results taken from NRQ^LUffiU. Using the 
equilibrium Green's function method in Sec. II, we de- 
rive an analytic formula for the the dependence of the 
self-energies in the Coulomb repulsion in Sec. III. We 
see that as the Coulomb repulsion increases the subgap 
states come closer to the Fermi surface and at a critical 
repulsion U c they overlap it. At stronger repulsion the 
states move away from the Fermi surface toward the gap. 
Bauer et.al. citebauer employed the NRG method, which 
is a reliable approach to analyze low temperature physics 
of spectral functions. The border between singlet and 
doublet ground states (i.e., Tk/A rj 0.3) was estimated 
by the crossover of bound states from the Fermi surface. 
We reproduce this result and present an analytical for- 
mula for this border. Within the transparent transmis- 
sion regime, the supercurrent is worked out in Sec. IV 



and it is shown it changes its flow direction in strong in- 
teractions. In the final section, the current noise due to 
the charge blockade in the junction is worked out and 
the microresonators are detected for a different Kondo 
temprature. 



II. SINGLE IMPURITY INSIDE THE 
JOSEPHSON JUNCTION: THE GENERAL 
FRAMEWORK 

A. Hamiltonian 

Consider two s-wave superconducting reservoirs j = 
L, R, different only in their phases, described by 
the mean-field Hamiltonian W := J2k s e k n k s + 

J2 k Aexp(i^) 4 T ci fc i + h.c, where n{ s := 4,14,8 is 
the number of particles with momentum k and spin s 
at the reservoir j. One magnetic impurity is inside the 
junction with the Hamiltonian 'H lmp := edn s —U n^n^, 
where n s := dld s denotes the number of particles with 
spin s. and the Coulomb repulsion U between electrons. 
The coupling between the impurity and the reservoirs is 
defined by H 3 T := tj^k s (4 s d s +h..c.) where t is the hop- 
ping matrix element. The total Hamiltonian becomes 

2 

n := H imp + (W + ?4) (!) 

i=i 

We assume the chemical potential is the constant Fermi 
energy (set as a reference energy), i.e., there is no bias 
voltage. We also consider the two reservoirs to be identi- 
cal on all of their parameters except of their phases, their 
phase difference is <f>x — —4> 2 = <j>/2. We further assume 
particle-hole symmetry = —U/2. We also define the 
total coupling strength to the j-th reservoir Tj :— np\tj | 2 , 
where p is the density of states. 



B. Green's function 

We are using Green's functions in Nambu space^, here 
and henceforth marked by a hat " " ". The Green's 
function corresponding to the Hamiltonian (JlJ in the 
non-interacting limit at the impurity site is knowrJ^l to 
be Gq(w) = (a) — €d<r z + o^T) -1 , where uj :— ui(l + 
J2i^i/E(u>,Ai)) and f := ^ TjAe^ /Eju, AQ. The 
function E(u>, A) — —i sgn(o;)v / w 2 — A 2 above the gap is 
the quasiparticle energy, and for the (imaginary) frequen- 
cies below the gap as E(uj, A) = \/A 2 — uj 2 . With the 
choice of identical superconducting reservoirs one can see 
Co = u(l + 2T/E{u,A)) im&T = 2TAcos{<j>/2)/ E{lu, Ai). 

The Dyson-Gorkov equation formally gives the inter- 
acting Green's function in terms of the diagonal and off- 
diagonal self-energies E and Ea 
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Figure 2. (Color online) Energy diagrams illustrating Cooper pair transport through a magnetic impurity, (a) Weak Coulomb- 
interaction limit where transport occurs through a screened spinless level; denoted by screened arrow at the bound-state Et- 
The Cooper pair parity is conserved, (b) Strong- Coulomb-interaction limit where transport flips the impurity spin. The parity 
of the Cooper pair is reversed, which results in a reversal supercurrent flow. 



A 1 (S>-e d - S(cj) f-E A (0) \ 

d F(u,<t>) \ f*-E^(0) w + e rf + E*(- W ) J 

(2) 

where F(ui) is the determinant of the matrix. 

C. Going beyond mean field 

As one analytical approach, we are choosing second- 
order perturbation theory in U /2nT. It will provide 
analytical expressions for the subgap bound-state ener- 
gies and their spectral weights, which are otherwise pre- 
cisely known only numerically through NRG results. The 
lowest order mean field theory, on the other hand, has 
the known restrictions already discussed. We will see 
that from finding the second-order self-energy matrix in 
Nambu space (and in the particle-hole symmetry), we 
can extract analytical formulas for the bound states and 
their spectral weights in the non-magnetic phase. We 
noticed these results are in agreement with NRG results 
and thus reliable for our purpose to study noise in the 
domain of validity of this perturbative approach. 

The significant improvement from including second- 
order perturbative theory at least in the non-magnetic 
phase relative to mean field theory can be tracked down 
to a set of reasons. The self-energy in the particle- hole- 
symmetric case has (diagonal) particle-particle as well 
as (off-diagonal) particle-hole matrix elements in Nambu 
space. In the mean field approximation the Coulomb re- 
pulsion modifies the diagonal self-energy by shifting the 
bare electron and hole energies and the off-diagonal self- 
energy by inducing an order parameter into particle-hole 
pairs. This self-energy leads to poles of the Green's func- 
tion, i.e., induces bound states below the superconduct- 



ing gap. The subgap en ergies were found analytically^ 
as well as self-consistentlySSED The second-order pertur- 
bation adds a term that depends on U 2 , T, and ui into 
the diagonal self-energies and a term that depends on 
U 2 , r, and <j> into off-diagonal self-energy. Some of the 
improvements these terms cause are listed below. 

1. The diagonal mean field self-energy is frequency- 
independent P^l This leads to a vanishing 
Kondo temperature Tk, i.e., it is insensitive 
to interaction — see next section. Lacking an 
interaction-controlled Kondo temperature the 
singlet-doublet transition, which is supposed to be 
governed by Tr-/A, becomes uncontrollable. The 
second-order self-energy improves this by providing 
a frequency-dependent self-energy and therefore 
an interaction-sensitive Kondo temperature — see 
section IIIII 

2. Accurate NRG results show that, as the repulsion 
between electrons increases, the bound states come 
closer to the Fermi surface, and after overlapping it 
at a critical interaction U c they slowly move away 
from the Fermi surface. This behavior is captured 
in the self-consistent mean field solution^, how- 
ever, the criticality appears to be independent of 
the superconducting gap A. This contradicts the 
precise NRGP3 result as U c must decrease for su- 
perconductors with larger gaps. The second-order 
perturbation improves this into an NRG-like A- 
dependence of the criticality. This dependence can 
be seen in Eq. (JsJ below. 

3. The phase boundary between the singlet and dou- 
blet ground states in the strong coupling regime 
was worked out by NR G methods and estimated 
by T K /A w Q.3) 59|74 | 75 | This boundary can also be 
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approximated by the criticality line IfJ^ However, 
in the mean field, the Kondo temperature is not 
tuneable. In the second-order the border Tk is 
recovered as a an inversely scaling interaction for 
which this border can be well-approximated by the 
criticality line — see Fig. 

4. The self-consistent mean field solution describes su- 
percurrent reversal in the strong Coulomb repulsion 
regime, a — tt transition. This transition turns out 
not to be sensitive to the relevant physical param- 
eters, in contrast to more accurate s olutions sug- 
gesting that it should be controllable! 59 1 60 1 63 ^ This 
sign change is recovered in a controllable way in 
the second-order with a modified current-phase re- 
lation due to the cp dependence of the bound-state 
energies appearing in that order. 

Let us emphasize our solution is analytical and re- 
stricted to the non magnetic phase. Also we are not 
aware of any self-consistent analytical solution in the en- 
tire phase diagram. 



III. SECOND-ORDER PERTURBATION 
THEORY 

A. Self-energy 

The second-order perturbation of the Coulomb inter- 
action in the Nambu space can be computed by Feynman 
diagrams using normal and anomalous propagators. The 
diagonal self energy in the first order is known to be^ 
£W = U/2, which shifts the bare energy of impurity 
into id = e + U/2. Usually the particle hole symmetry 
is imposed onto the problem to set £d = 0. The off- 
diagonal self energy in the first order is = — A^, 
where is the induced superconducting order parame- 
ter at the impurity site. By definition, A^ := U(d^d^) = 
^ J °° duj Tr (a x G d (uj)), which contains the off-diagonal 
element of the impurity Green's function. 

As discussed in the introduction the low energy physics 
is sensitive to the Kondo temperature Tk which essen- 
tially measures the width of the subgap many-body res- 
onance. The Kondo temperature in perturbation theory 



is defined as Tjc 



l-dz/duj a t w = 0.^ In the sym- 
metric Anderson model this temperature vanishes in the 
mean field limit. Therefore, any phase transition in the 
junction that is governed by the ratio of Tk/ A, if even 
occurs, is unphysical. 

In the second-order , the self-energy is defined as 



£([/) 



a(U) and S A := P(U,<f>) A, where a(U) 



a{U)uj, and (3(U,cj)) = + ^ (U, <j>) for pW(U,(f>) = 
b(U) cos(</>/2). These formulas are extracted from Feyn- 
man diagrams and details are given in Appendix A. Con- 
sidering this order of perturbation theory the Kondo tem- 
perature not only does not vanish, but also it scales in- 
versely with U, which although fails to provide the exact 



exponential decay, yet supports decaying with increasing 
Coulomb repulsion. 

Note that in the particle-hole symmetric case the 
Friedel sum rule remains satisfiedP3 . However, beyond 
this symmetry a modification is needed on the self-energy 
in order to satisfy the Friedel sum rule. A typical modifi- 
cation was pr opose d by Martin-Rodero and Levy-Yeyati 
and coauthor d ' " ' that maintain the sum rule. A brief 
overview of this modification can be found in the Ap- 
pendix B. Nevertheless, because in this paper we consider 
the particle-hole symmetry for granted, the Friedel sum 
rule is satisfied and no modification is required. 

The second-order Feynman diagrams read 



t) (2^) 



(3) 
(4) 



The diagonal self-energy was first reported in 1975 in 
normal metals by Yoshida and Yamada^^ in the An- 
derson model in the normal conducting state. As dis- 
cussed earlier in Ref. ITU and mentioned in Eq. ( Al I the 



normal propagator in the small-gap limit can be assumed 
not to be modified by superconductivity. The supercon- 
ductivity effect enters the off-diagonal propagator. Hav- 
ing written the effective Hamiltonian as the summation 
of non-interacting part and the self-energy interacting 
part, the perturbative energy ratio of diagonal Hamil- 
tonian is almost ~ Eb(U)/8A in the unit of (U/ttT) 2 . 
This indicates that the perturbatively relevant limit is 
U/ttT < 2.8WAjEf,, thus the closer a bound-state is 
to the Fermi surface the more accurate is the pertur- 
bation theory on diagonal interactions. However, for off- 
diagonal interactions the perturbation limit turns out to 
be U/ttT < 2.5. It is important to emphasize that the re- 
sponsible physics can only be inferred within the limit of 
min{2.5, 2.8y/ A/Eb}. Note that in this section, because 
we want to compare our results with NRG results, we 
go a bit beyond the perturbation limit and interestingly 
reproduce all NRG physics. Nevertheless we will remain 
in the perturbative limits for the rest of this paper where 
supercurrent reversal and its decoherence are discussed. 

Note that current in a not very long weak link van- 
ishes at the phase d> = Tt which is also confirmed in 



experiments This would not be achieved by the self- 
energy found in Ref. |76| . but by the one found in this 
paper. 



B. Bound-states 

The magnetic moment at the impurity site is de- 
fined by the difference in occupation between the 
spin polarizations, m := (—U/2)((rif) — (nj,)) = 
~¥rr Io° ^ Tr(<7 z G d (w)). This contains the difference 
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between two diagonal elements of impurity Green's func- 
tion. This definition describes the spontaneous symme- 
try breaking responsible for distinguishing magnetic from 
non-magnetic states. Because of the possible difference 
between spin polarizations we define rj :— (e^ + U /2) — m 
as the (bare) energy of magnetic impurity. However, we 
restrict this study to the non-magnetic phase sector, thus 
m = everywhere. Due to the particle-hole symmetry 
the first order diagonal self-energy becomes U/2. There- 
fore in Eq. ([2]) the bare energy is summed up with 
the first order self-energy to be renamed r\. This shifted 
energy is now playing role as a new bare energy to be 
renormalized by the second-order self energy. 

The subgap states appear where the denominator of 
the Green's function in Eq. Q vanishes, i.e., F(±Eb) = 
0. The bound-state energies of these states are the solu- 
tions to 



(5) 

where E 2 := A 2 — E 2 . In the particle-hole symmetry and 
non-magnetic sector 77 = the bound-state equation is 
reduced to 



E bT A cos \ + {E b ~ a(U)E b ± (3 (U, 4>) A) ^ p = 0. 

(6) 

Let us briefly discuss the phase dependence of the in- 
duced superconductivity. We recall that by definition (3 
is Ad/ A to first order. Usually in non-self-consistent so- 
lutions, such as Matsumoto in Ref. (SUJ one 'chooses' the 
phase dependence of A^ such that a smooth phase depen- 
dence of the localized excited states is achieved. However, 
the exact phase dependence of A<j can be determined by 
solving self-consistency equations. In the high transmis- 
sion case that we study here, we assume that the phase 
of the induced gap at the impurity site is governed by the 
two reservoirs and therefore one can choose A^ to have 
a cos(0/2) dependence on the phase. This allows the 
P function to save this phase relation up to the second- 
order . On the other hand, the induced superconducting 
correlations at the impurity site are expected to drop 
through the repulsive Colulomb interaction therefore one 
can expe ct Ag to decrease with the increase of Coulomb 
repulsion^ 21 . 

From Eq. |6| one can easily deduce that the only al- 
lowed bound-state at (j> = tt is Eh = 0. Therefore, we 
consider the ansatz E b = Af(U, <j>) cos</>/2 as the bound- 
state in a high transmission channel. Substituting this 
ansatz into Eq. (|6| the function f(U, 4>) has to be the 



solution of 

(cos 2 ^(l-a) 2 ) / 4 ±^cos 2 ^(l- a )) f 

+ ^ + 6 W|-(l-a) 9 ) f 

± (-2 7 2 - 2b + 2ab) f + 7 2 - b 2 = 0. (7) 

where 7 := 217 A. 

Fig. (|3^) presents the sensitivity of bound states to 
the Coulomb interaction U, extracted from the solution 
of Eq. ([7]) . In the domain of second-order perturbation 
theory these states coincide with the numerical results 
taken recently from NRG method- 11 . Here we use those 
parameters used in Fig. (3) of Ref. [21] — with a remark 
that the total hybridization strength with both reservoirs 
+ Tr = 2r in our work is equivalent to the hybridiza- 
tion strength T in Ref. [51J where a single reservoir was 
considered. 

By the increase of the Coulomb repulsion the bound 
states first approach the Fermi surface and, at the critical 
repulsion, overlap with the Fermi surface. Moreover, one 
can see in Fig. (|3^) the larger the gap is, the smaller 
the critical Coulomb repulsion becomes. Note that this 
A-dependence of the critical repulsion can be seen only 
after the second-order is involved in the self-energy. This 
is the reason why this phenomenon cannot be seen in 
the first-order studies, neither in analytical nor in self- 
consistent first-order solutions.^! 

The most relevant regime for the parametrization is the 
large gap regime. We can study this regime by param- 
eterizing our bound states in terms of A/2r and A/U. 
The resulting bound states are plotted in Fig. ([3|. 

The critical Coulomb repulsion can be found by setting 
bound states to zero in Eq. ^. This gives rise to the 
analytical formula for the critical repulsion 

/3(C/ c ,^)= 7 cos^. (8) 

where 7 has been defined below Eq. (m). By substituting 
the off-diagonal self-energy into Eq. p) one can find the 
following equation whose solution gives rise to the critical 
Coulomb repulsion in terms of A and F. 

(y-l) (^Acos^ + C/.A^rcos^. (9) 

As one can see in Eq. |9]) for fixed T and increasing A 
the critical repulsion U c decreases as A -1 / 2 . This behav- 
ior is in qualitative agreement with NRG results in that 
U c grows with falling A at roughly the same rate, albeit 
data are not sufficient to compare with the power law in 
detaiP. 

As U increases and crosses U c the spin orientation sym- 
metry is broken and a magnetic moment forms, which 
indicates the magnetic impurity becomes less screened. 
This will introduce a nonzero 77 in Eq. (pi. The bound 
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Figure 3. (Color online) Bound states at (a) zero phase limit cj> — due to the localized magnetic impurity inside a Josephson 
junction between two superconducting reservoirs, where A on the solid line is 0.06, dashed line 0.01, and dotted line 0.001; (b) 
in variable phase < <j> < tt for different M^, where A = 0.06. In both plots 7iT = 0.01. (c) The bound-state in terms of A/r 
and A/U for the zero-phase case (<f> = 0). 



states in the magnetic phase can only be found through 
the self-consistent analysis of the problem as the mag- 
netic moment (say rj) is determined at the same time 
as the bound states are. However in the non-magnetic 
phase we noticed because the bound states must be zero 
at <f> — (odd number) x it the bound states are of the 
form f(U, 4>) cos(0/2), or in other words f(U, 4>) (1 — 
r sin 2 4>/2) 1 / 2 with r = 1 for a "ballistic" channel. Beyond 
the non-magnetic phase one can expect the transmission 
rate r to decrease with the rise of the magnetic moment 
in the junction and to depend on the spin polarization at 
the impurity site. It is natural to seek a solution in the 
magnetic phase, however, due to the inherited restriction 



of solutions in perturbative approaches in strong inter- 
actions we do not consider non-ballistic transmission in 
this paper. 

C. Spectral weights 

The weight of the bound states is determined by eval- 
uating the residue of impurity Green's function at the 
bound states. The normal (diagonal) and anomalous (off- 
diagonal) weights are denoted by W and Wa ■ Having de- 
fined z(U)~ l = 1 — a(U), one can find the weights from 
Eq. ^ as 



E_ +z± vzE_ 

W(±E b , U,4>) = ^ ^ 2TEb - — — (10) 



1 + Wz + W + W z P(<t>, U) cos f + ^zsm' f 

w A{ ±E b , u, a = ±^ E2 

2Eb 1 + & + 5 + U) cos t + -n 2 t 



where E = ■\/A 2 



E%(U, 4>). Note that in Eqs. (|10| and 



(11 1 E and E b depend on both U and (j), and z depends 



on U. 

In the particle-hole symmetric case one can see that 
the normal weight is an even function of energy and the 
anomalous weight is an odd function, i.e., W(Eb, U, (f>) — 
W(-E b ,U,ct>) and W A {E b ,U,<t>) = -W A (-E b ,U,&. 
The asymmetric case needs a modification of the defi- 



Bl) 



nite self-energies to become a and /?, (e.g. see Eq. 
in Appendix B]. 

The diagonal and off-diagonal weights associated with 
those bound states whose energies were already given in 



Fig. (|3ja) are indicated in Fig. Qa) and[4jb), respec- 
tively. The diagonal spectral weight for A = 0.06 and 
0.01 are in good agreement with the corresponding NRG 
results however, the criticality of A = 0.001 C ctSG j ctS 
indicated by NRG, is ~ 57TT which is beyond the validity 
of perturbation theory. 

The off-diagonal weight W A changes its sign at U c . 
Due to the proportionality of U c and A -1 / 2 (via Eq. 
one can expect that in larger gaps this change of sign is 
shifted to weaker interactions. In next section we will see 
that the supercurrent flow inside a junction depends on 
the off-diagonal weight and therefore one expects that 
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with decreasing Kondo temperature to below 0.3A the 
supercurrent reversal occurs. 

Moreover, as the gap becomes smaller the maximum 
weight is shifted to higher Coulomb repulsions. Namely 
when the ratio of 7 is about 1 the maximum weight is at 
zero interaction and decreases monotonically in stronger 
repulsion; however for larger 7 the maximum is pushed 
to somewhere between zero and the critical repulsion. 

The competition between superconducting gap A and 
Kondo temperature Tx becomes more clear in a phase 
diagram with the background color being the off-diagonal 
weight of the bound states. For the special case of zero 
phase this is depicted in Fig. ([5]). 

In this phase diagram the dotted line describes the 
phase boundary of Eq. ([9]|, defined as the line where 
the bound states overlap the Fermi surface. This coin- 
cides with the boundary where the off-diagonal weight 
changes its sign. The off-diagonal weight as discussed 



below Eq. (Ill should become negated as soon as the 



bound-state crosses the Fermi surface therefore it must 
either smoothly or suddenly change its sign about this 
point. Let us consider a point very close to the Fermi 
surface = e where f3 ^ 7COs(</>/2). Expanding the 
bound states from Eq. (fsl) about this point one can find 
e/A = q=(7Cos(0/2)— /3)/(7+l— a). Substituting this en- 
ergy into the off-diagonal weight one finds the numerator 
is finite valued in the vicinity of the Fermi surface while 
the denominator because of the presence of Ej, becomes 
very tiny. This indicates that the off-diagonal weight 
should jump discontinuously from positive to negative 
when crossing the Fermi surface. 

This can be seen from Eq. (Ill where Eb = both 



the numerator and the denominator approach to zero. 
The off-diagonal weights, the pair amplitude, measure 
the degree of induced superconductivity. They connect 
to Ginzburg-Landau theory by taking the limit of inter- 
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Figure 5. (Color online) Phase diagram for singlet and dou- 
blet ground-states. The critical Coulomb repulsion is denoted 
by U c to be compared with the doublet-singlet boundary at 
Tx/A = 0.3. The dashed vertical line is the doublet-singlet 
boundary at large A limit. The dash-dot line indicates the 
border between singlet and doublet states for very large gap 
superconductors. The off-diagonal weight Wa/A is plotted 
in background colour. 



n,j,H]/th. By substituting the Hamiltonian 



of Eq. |1| and denning Ae 1 ^ ~ Y.k^k.^-k.^ the 



current operator depends on impurity-reservoir Green's 



+1 ° function G dj , i.e., J = (2et/h)hnf™ du)Tr(G J > (u>) - 



G>(<jj))- Due to the Dyson equations the impurity- 
reservoir Green's function depends explicitly on the im- 
purity (Q d ) and the reservoir (g l ) Green's functions, i.e., 
namely G r j /a (u) = tg J r/a (uj)a z G d / a (uj). The Schwinger- 
Keldysh Green's functions are determined by the re- 
tarded and advanced Green's functions G>{uS) = (1 — 

f{u))){Gr{u)~Ga{u)) and G<(u) = -f(u)(Gr(u)-Ga(u)) 

where, in thermal equilibrium, / is the Fermi function 
at temperature T. Superconducting Green's function 
is known to be (g J )n = (g J )22 = irpuj/ E(uj, Aj) and 



12 



21 



npAje l Vi /E(u,Aj), where p is the 



density of states at the reservoirs. One can easily show 
at finite temperature that the current density between 
the impurity and the j-th reservoir becomes 



h 



[g 3 (u) 21 G d (u) 12 ] f(u)du. (12) 



action constant gees to zero, [i.e., >Fe 1 ^ = ffBCs(^t^i) m 
a superconductor]. Sign change in the right side can be 
interpreted as the phase shift by it in the pair amplitude. 
This phase lapse is traditionally called a — n transition. 
In the next section we will discuss this — ir transition 
through the supercurrent reversal flow. The dash-dotted 
line is the boundary at which the ground state makes 
a singlet-doublet transition in the large- A limit, (i.e., 
U/T = 2). Singlet ground states [i.e., (|) + | \] r ))/y/2\ 
are on the left side of this border and the right side 
contains doublet groun states (i.e., {| f), | |)})]S11Z6J The 
dashed line gives the transition governed by the Kondo 
temperature T K = 0T 82UV8r/7rE/exp(-7rL//8r). Nu- 
merical analysifPEHzH indicates that T K = 0.3A is ap- 
proximately the phase boundary between magnetic and 
non-magnetic phases. This line is almost in the nar- 
row region between the U c border and the maximum off- 
diagonal weight. 

One can see in Fig. ^ that the critical Coulomb repul- 
sion line approaches the phase boundary line Tk /A = 0.3 
in the large U. This indicates that in large U limit the 
crossover of bound states occurs in the vicinity of singlet- 
doublet transition 



IV. SUPERCURRENT 

A. Definitions 

The Josephson current in a junction is derived from 
charge conservation fid + Jl — J_r = where rid = til + nrt 



Without an impurity and identical reservoirs this cur- 
rent turns out to be sinusoidal with respect to the phase 
difference <f>. However, in the presence of a magnetic im- 
purity this sinusoidal current-phase relation is deformed 
due to the presence of subgap states. Note that within 
the Hamiltonian model of Eq. (JTJ the Josephson current 
is conserved: (J L ) = —(J R ). 

The impurity Green's function describes the discrete 
bound states in the subgap region as well as the contin- 
uum spectrum above the gap. The bound states appear 
as poles of the Green's function. Accordingly, the Joseph- 
son current has two contributions: 1) from the contin- 
uum, J J C , which appears by tunneling into supergap states 
at the impurity site, and 2) from bound-state current 
that appears due to the tunneling of electrons into sub- 
gap states. Some collective properties of both have been 
studied using NRG approach by Choief al. 2004.1221 More 
recently, Karrasch Oguri and Meden studied the total 
supercurrent by NRG and compared with the functional 
renormalization group (fRG) methods^. As a result the 
total supercurrent manifests a rather complicated de- 
pendence. With the increase of U the current starts to 
be suppressed, until it starts to flow backward near the 
critical U c . 

Mean field theory can also predict the current reversal, 
however it lacks the Kondo temperature as discussed in 
Sec. I. Let us now study the current flow in the second- 
order . We will see that the current flow can exhibit a 
controllable flow in the reverse direction in the strong- 
interaction limit. 
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B. Subgap contribution to the supercurrent 

tunneling of Cooper pairs inside a junction between 
two superconducting reservoirs gives rise to the free 
energy of the junction U(<fi) = — Ej cos 4>, where 
Ej is the Josephson energy. A nonzero value of 
<j) will cause a supercurrent I sc = (2e/ti)dU/d4> = 
(2eEj/h)sui(j) across the junction. In a tunnel junc- 
tion without traps or time-reversal symmetry-breaking 
impurities, the Ambegaokar-Baratoff formula^ predicts 
Ej = (/iA/4e 2 i?)tanh(A/2T), where R is the normal 
state resistance of the junction and T < T c . At nonzero 
temperature due to the thermal activation of quasipar- 
ticles there is also the normal current in addition to the 
supercurrent. However, at small temperature T <C 2Ej 
the effective resistivity of the tunnel junction becomes 
vanishingly small, thus the normal current becomes neg- 
ligible and the dominant current will be the supercurrent, 



i.e., I <C isc- This supercurrent is made of two parts: 
the continuous state and the subgap state parts: I sc — 
^subgap + icont. From the Ambegaokar and Baratoff for- 
mula one can expect that for temperatures much smaller 
than one half of the superconducting gap T <C A/2, the 
continuous part becomes negligible, thus the subgap part 
becomes dominant. Here, we work out the subgap part 
of the supercurrent, which we expect to make the most 
crucial contribution. 

The bound-state contribution sum up over all bound 
states. At zero temperature there is only the contribu- 
tion of subgap states below the Fermi surface, there- 
fore by substituting the off-diagonal spectral weight, 
instead of full Green's function, in Eq. ( 12 1 one 



finds the subgap contribution to the current as {Ij) = 



{l6eTA/h)W A sinO/2)/£, which by the use of Eq. ([TTl 
turns into 



L 



4erA 2 



sin 



subgap 



<t>-m<t>) 



fsinf 



E b E jjZzP(4>,U) cosf 



A 2 
E 2 



E 
lYz 



tanh 



\E h 



2k B T 



(13) 



where E and Eb depend on both U and <f>, and z depends 
on U. 



The subgap supercurrent dependence on </> at T — 0.1 A 
is shown in Fig. (|6^) for different U . One can see the 
suppression by stronger repulsion and the change of sign 
at U c . Fig. (J^Jd) shows the current reversal as a func- 
tion of repulsion for different values of the phase at zero 
temperature. One can see that the extreme currents flow 
about the phase 7r. As indicated in Fig. ([6^), by the 
temperature rise, the maximum supercurrent, whose ex- 
act value depends on the U and temperature, is shifted in 
the phase into between zero and n. By increasing repul- 
sion the extrema are shifted from nearby n toward tt/2 
and after the crossover it is pushed back again towa rd 7r. 
This behavior has also be seen in NRG studiePHSl. 



Note that the functional renormalization group in Ref. 
[63] shows that the supergap states may contribute to 
the (j> dependence of critical repulsion U c . As a result 
the supercurrent reversal depend on both U and (f>. In 
other words, the criticality occurs for a particular re- 
pulsion at a particular phase. Let us emphasize that 
here we restricted our analysis to the subgap states con- 
tributions and therefore our analysis lacks this peculiar 
phase dependence. Nevertheless, what practically can be 
measured is not the exact supercurrent at a particular 
phase but instead the average of supercurrent at differ- 
ent phases. 



V. CURRENT BLOCKING DECOHERENCE 

Spectroscopy experiments on phase and flux qubits 
spectroscopically reveal the presence of a few microwave 
resonators! 68 ! 84 ! These microresonators behave as spuri- 
ous two-level systems inside the junction, whose cou- 
pling to the qubit produces reduced measurement fi- 
delity and decoherenceP^. They can be intrinsically very 
coherent JHsHHU Recently, the trapping scenario for the 
dipole-dipole interaction between these states and the 
conduction electrons was studied by Constantin et 
Using a different method, deSousa et.al. citedeSousa con- 
sidered the noise spectrum due to the occupation of a sin- 
gle localized magnetic impurity hybridized with electrons 
in a superconductor in the non-interacting regime. Their 
study shows the presence of some resonances in low tem- 
perature that could serve as candidates for the observed 
microresonators. However, the Kondo temperature in the 
non-interacting limit vanishes and their study does not 
include the competition between the Kondo temperature 
and the superconducting gap. On the other hand Faoro 
and Ioffe found that the underlying 1/f noise spectrum 
behaves correctly with the temperature [i.e., S(ui) ~ T 2 ] 
only if the subgap bound states can be formed in the 
vicinity of the Fermi surface, which is not the case in 
the non-interacting Andreev bound states. As it was 
explained earlier in the paper the subgap bound states 
come closer to the Fermi surface if the Coulomb repul- 
sion is taken into account. Therefore, Faoro and Ioffe 
concluded that the underlying 1/f noise must originate 
in interacting electrons. Here, we study how the spec- 
trum of additional microresonators may be modified in 
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Figure 6. (Color online) Supercurrent variation with respect to (a) phase difference in different Coulomb repulsions at kT = 
0.1A, (b) the Coulomb repulsion U at different phases between and 7r at zero temperature; the reservoir-junction hybridization 
rate tiT = 0.1 and the gap A = 0.01. I c = 



the presence of interaction.^^ We will see the range of 
second-order perturbation theory is sufficient to observe 
the emergence of novel effects in the noise spectrum as 
the consequence of competition between superconducting 
gap and the Kondo temperature. We should note that a 
scenario based on local dipoles instead of traps has also 
been brought forward^. 

In general the current noise spectrum is defined as 
S = J™ dt([5I(t),5I(Q)]+). One can consider a con- 
stant hybridization coupling between the reservoirs and 
the junction and substitute the critical-current into this 
definition. The outcome from each magnetic impurity is 
known to be a Lorentzian noise spectrum, and a combi- 
nation of several such spectra leads to 1/ f noise. 



For the purpose of understanding the nature of a few 
microwave resonances on top of this noise spectrum one 
can assume that the magnetic impurity partially block 
conduction whenever it captures electrons from a reser- 
voir. The model proposed by deSousa et.al. citedeSousa 
assumes that the channel average matrix element for elec- 
tron tunneling from one lead to the other depends on 
n according to t ~ + t^n. From Eq. (12) this 
can vary the supercurrent as I c — J c (l + 
This modifies the current noise spectrum by the varia- 
tion of number of electrons at the impurity site n such 
that StotaiM = Sj((j) + {5I c ) 2 S n (uj). We shall study 
the contribution of the second term in the presence of 
Coulomb repulsion. This noise for u > is 



/CO 
de[r CTCT '; CT v(e ; e - w;e - cj, e)A tTl<J >(e)A <T><7 >{e - u) 

(7,(7' -°° 

-r'crcr'^'o-O, e-oj;e-u;, e)B a ^ a {e)B a > , a (e - u)] [1 - /(e)]/(e - uj) 



(14) 



where r <7CT / ;CT / .(ei, e 2 ; e 3 , £4) and T'„i a . a<T > (ei, e 2 ; e 3 , £4) are 
the diagonal and off-diagonal vertex corrections, respec- 
tively. The diagonal (off-diagonal) density of states A (B) 
consists of the normal (anomalous) weight W (Wa) at 
the subgap states and the continuous function A (B) for 
the supergap states. An analytical formula for the con- 
tinuous part of the density of states is, to the best o f our 
knowledge, unknown. However, numerical analysis^ 43 ! 51 ! 



indicate that with the increase of U the continuous den- 
sity of states above the gap is suppressed about the gap 
while two smooth maxima appear in high frequencies. 

The noise spectrum has three major contributors in 
three frequency ranges (1) the subgap-subgap spectrum 
at low energy, (2) the subgap-continuous spectrum at 
moderate to high frequencies, and (3) the continuous- 
continuous spectrum in high frequencies. As qubits are 
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Figure 7. Vertex correction in the second-order 



operated at frequencies much below the gap, the first 
contribution is most relevant. For the purpose of study- 
ing the microresonators that appear on top of the 1/f 
noise spectrum we compare the weights and frequency 
of peaks in the interacting noise power spectrum with 
that of non-interacting model 42 . The vertex correction 
diagrams are shown in Fig. The diagrams (a)-(c) 

represents the normal and diagrams (d)-(f) the anoma- 
lous contributions into the second-order interaction of 



vertex correction. Analytical calculation of these dia- 
grams shows that the diagrams (b) and (c) are exactly 
canceled out and so do the diagrams (e) and (f). More- 
over, the diagram (d) has a tiny contribution in the limit 
of A < r which is negligible. One can calculate the ver- 
tex correction functions T » 1 + (r 2 7r 2 /8)([//7rr) 2 and 
V « 1. The noise spectrum at low frequencies can be 
separated in the following form: 



S n (u)/2h= 5{u-2E b ) 2W 2 T[l - f(E b )]f(-E b ) 

+ 9(lo - (A - E b )) {W[A{-E b -u) + A(E b + u))T - 2W A B(-E b - u)}[l - f(-E b )]f(-E b - u) 
+ 9(lj - (A + E b )) {W[A(E b -u)+ A{-E b + u))T - 2W A B(E b - u)}[l - f(E b )]f{E b - u). (15) 



Note that with the increase of U the subgap bound 
states approach the Fermi surface. Within our focus of 
low frequencies two noise spectra are depicted in Fig. 

Fig. ^ presents the charge blockade noise spectrum 
for various Coulomb repulsions. With the increase of 
Coulomb repulsion the spectral noise of subgap-subgap 
transition is affected in two ways: 1) the resonance fre- 
quency is lowered until it becomes zero at the critical 
repulsion. Stronger repulsions than the critical one are 
expected to push back the resonance into higher frequen- 
cies. 2) The noise spectral weight of the sub-gap reso- 
nance by the increase of repulsion decreases (except at 
small gaps where the noise may first increase and after 
achieving a maximum at U < U c it suppress onward). 
As a typical example one can see that with the increase 
of Coulomb repulsion from zero to U = 0.2 the weight 
of the resonance is suppressed by 60%. We recall that 
the normal weight, shown previously in Fig. Qa), in a 
large-gap superconductor — i.e., A ~ T monotonically 
decreases with the increase of U. However, by increas- 



ing U in smaller gap materials this weight may raise to 
a maximum at U max < U c and then becomes suppressed 
in stronger interactions. This can be seen for two typical 
parametrization of A/T = 0.94 and A/T = 0.28 in Fig. 
(J8^a) and Fig. (|8jb), respectively. 

The resonance noise is highly sensitive to the size of 
the superconducting gap, the hybridization strength, and 
the Coulomb repulsion between onsite electrons. The 
peculiar sensitivity of the spectral weight of this noise is 
plotted in Fig. (|9| for the two cases of <j> = and (f> = tt/2. 
The dotted line indicate the singlet /doublet boundary. 
The background color indicates S(oj)/A 2 with red the 
highest noise power that decreases into the lowest purple 
area through yellow, green and blue colours, respectively. 
In Fig. (|9| the left plot represents a zero-phase junction. 
The phase difference in the right plot is instead tt/2. By 
comparing these two one can see in the latter case that 
the power of resonance noise monotonically suppresses 
with the increase of interaction in every parametrization 

ofA/r. 
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Figure 8. (Color online) The charge blockade noise spectrum for various Coulomb repulsions as listed in the legend 




Figure 9. (Color online) The background colours indicate the spectral weight of the resonance noise as a function of ratio A/T 
and U/ttF for the phase difference = (left) and <j> = 7r/2 (right). Dotted lines are the singlet-doublet boundary. 



VI. CONCLUSION 

We presented an analytical formalism to study a mag- 
netic impurity residing at the Josephson junction par- 
tially blocking the only ballistic channel of the junction 
in the limit of moderate Coulomb interactions and small 
superconducting gap. We applied the second-order per- 
turbation theory in terms of the Coulomb repulsion and 
analytical formulate the subgap bound-state energy and 
their normal and anomalous weights. These results were 
shown to be in good agreement with recent results of 
NRG. We expressed the supercurrent ctS £l function of 
the Coulomb repulsion and noticed that in the strong- 
interaction limit the supercurrent may flow in the reverse 
direction compared to the noninteracting case. We finally 



studied the noise generated due to the partially blockade 
of conduction caused by the magnetic impurity. Depend- 
ing on the fabrication of the superconducting gap, the 
hybridization strength and the Coulomb repulsion the 
resonant noise can be heavily suppressed from the super- 
conducting material. 

Note that we do not have control on the magnetic im- 
purity site. For instance the level position e e i, measured 
from the Fermi energy of the two leads is not a parameter 
that is tuneable by an external gate voltage, as it is in 
the case of quantum dot. The interaction U is not a con- 
trollable parameter either. As a consequence we cannot 
think of the bound states as a qubit. These states play 
the role of spurious qubits which could strongly couple 
to the desired qubit. However the microwave resonances 
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they generate for the critical-current are highly sensitive 
to the interaction. 

A quantum dot when coupled to two superconductors 
may behave as a magnetic impurity in the junction. One 
of the advantages of using a quantum dot is that the 
Coulomb repulsion is fully controllable in a quantum 
dot, thus one can test the suppression of the resonant 
noise spectrum. In this paper, we studied only the non- 
magnetic solution. Extending this study to cover the en- 
tire phase diagram, one needs to set up a self-consistent 
solution. 

The consequences of this physics of superconducting 
qubits can be elucidated as follows. It is established that 
junctions contain trap states like ours. Similar to the 
noninteracting case^l, depending on the trap energy, the 
junction resonators 68 can be anywhere within the super- 
conducting gap, hence most of them won't ever resonate 
with the qubit whose energy splitting covers a band that 
is deep inside the gap and comparably narrow. Our work 
shows that Coulomb repulsion moves these resonators to 
lower frequencies hence toward the qubit frequency, while 
simultaneously reducing their spectral weight. With suit- 
able distributions of both e ( i and U, we might expect 1//- 
type noise spectra in addition to the resonator noise. We 
also see that the 7r-junction mechanism, albeit intriguing, 
does not give a marked contribution to the noise. 

A quantum dot as an artificial impurity in between su- 
perconducting reservoirs can be one option to test our 
predictions. A quantum dot with an occupied electron 
breaks Cooper pairs to be screened and produce a Kondo 
singlet state. A variety of Kondo temperatures f or the 
states with odd number of electron are accessibl e 1 55 ! 90 !, 
where negative critical-current is expected. Having ac- 
cess to both distinguishing the quantum dot states and 
their associated Kondo temperature (i.e., the effective 
Coulomb interactions) the phase decoherence of an ini- 
tial state can reveal the underlying critical-current noise, 
which can be compared at different Kondo temperatures. 
The maximum decoherence is expected at a narrow zone 
of parameters. One can also think of applying a weak 
in-plane magnetic field to a phase qubit, or coupling it to 
a cavity in order to tune the Coulomb interaction in the 
junction and measuring the difference in the spectral den- 
sity of the resonance peaks in the excitation frequency. 



The dominant behavior of the normal Green's function 
in the small gap limits can be approximated by a normal 
metal: 



1 



D(oj) 



w 1 + 



r 



u + sgn{uj)Y 



(Al) 



where D(u) is the same as F(uj) denoted in the denom- 
inator of Eq. Q for the noninteracting case (U = 0). 
Note that in noninteracting particle-hole symmetry = 
0. 

The off-diagonal Green's function vanishes when the 
gap is set to zero; however at small energies there is a 
finite contribution that can be provided by taking the 
limit of uj — > before considering the small-gap limit. 
One can find this function as follows: 



Arcos(>/2) 



E(u) 



AT cos f 

w(w + r) 5 



(A2) 



Considering the Pauli exclusion principle the Coulomb 
repulsion takes effect between two electrons with opposite 
spins. In the limit of A <C T the following diagrams take 
part into the self-energy matrix: 



1. Diagram (a) 
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Appendix A: SELF-ENERGY 

Let us approximate the Green's function of a supercon- 
ductor with a small energy gap with simpler functions. 



From Eq.(2.8) in Yamada and Yoshida's paper^SHH] 
the send order term of G(t, t') is given by 



-U 2 {-\f 



{G(t,1)G(1,2)G(t 2 ,t') 



+ G(r, 2)G(2, 1)G(1, r')} G(l, 2)G(2, 1), 

where the term within { } is for up spin and 
G(l, 2)G(2, 1) is for down spin. Let is calculate the first 
term in the Fourier space. Using Feynman diagram laws 
one can write 
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As the states are all above and below the superconduct- 
ing energy gap, we consider the contribution of in-gap 
states separate from the integrals and sum the weight of 
these states with the integral. The following phase space 
diagram | A 1| shows the domain of integration that is 
restricted to the grey regions on top of the red lines (An- 
dreev bound states). The energy w b is the bare bound- 
state energy that can be reached by setting U = in Eq. 
(10 I. In the second line the second parenthesis is writ- 



ten considering negative y and the remaining lines are for 
positive y. By the use of the condition A <§; T the above 
integration is simplified into 




To simplify the analysis we considered the small 
energy-gap limit -c 1) and approximate A in the 
domain of integration to zero. Let us first do the contin- 
uous part. After a change of variables from y to — y in 
the first integral and with the change of variable in all 



integrals into s 
result is simplified to: 



M(lM+gD 



r 2 



for positive values of y, the 




Figure 10. The integration domain space. 



For small frequency uj we expand the integrands by the 
use of ~ 1 =F x. This gives rise to: 
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The first integral results into \ {^f U - ^ ) to. The 

second integral gives w-linear term as for small frequency 
lo. To see this we first expand the integrand and sim- 
plify it, then we take the integration that in the gen- 
eral case becomes 2 ln(l + s) j\J\ + s + 4/ \/l +~s — ln(l + 
y/l + s) ln(l + s) - 2dilog(l + y/1 + s) - 2dilogy / l + s . If 
we expand this for small s and calculate it in the integral 
domain from to ^ the result is proportional to s and 

the corresponding self-energy becomes | (^) 2 o;. The 
total contribution of the continuum into the self-energy 



becomes 



*) 2 ( 



3- 



■. The full self-energy con- 
tains both Andreev states on top of the continuum limit. 
One can show that the weight of bound states for the 
case of ojf, — A cos ^, via Eq. ( 10 ), W — 4 1 sin ^ | and 



Wa = ±W. Consequently in the limit of ^ the discrete 
contributions cancel out and the self-energy is left with 
its value at the continuum value: 



LO 



The two terms in Eq. ( A3 ) can be written for opposite 



spins therefore the 1/2 factor is doubled. 



Diagram b) 




SM = C/ 2 / g / ^G°(y)F°(LO + x + y)G°(x) 
In the leading this gives rise to 
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The third line is of second degree in A and we need 
only to compute the first two lines to represent the lead- 
ing order. One can compute after applying the approx- 
imation ln(l + x) sa x the second term in the limit of 

A < T, which gives rise to ^ J2 3 =+,- WaOV&Df anc * 
because Wa(lo) = -Wa(-w), the outcome is zero. In 
other words, there is no sin ^ contribution into the off- 



diagonal self-energy of Anderson model coupled to super- 
conductor leads. The central integrand in the first line is 
expanded with respect to lo. The term independent of lo 
is the only nonzero term for the limit of small lo because 
the term proportional to lo becomes an even function over 
the entire real numbers that is identical to zero. The only 
nonzero term is 



47T 2 2 



dy 
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dx- 



v - r J A y + r J 
i i 



(x + y){x + y + Ty x-T J A (x + y)(x + y + r) x + T 



dx- 



The result after taking the integration becomes: 



r 
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Appendix B: THE PARTICLE-HOLE 
ASYMMETRIC CASE 

To complete the discussion about the self-energy let 
us comment on the asymmetric case. In general one 



can break the particle-hole symmetry. In this case one 
should apply the modified self-energy computed from the 
interpolation between the small and large U irJ^ into the 
Dyson-Gorkov equation. This ends up with the following 
self-energy modification: 



E := 



l-9 2 (a 2 cu 2 - f3 2 A 2 ) 



8(a 2 oj 2 
(3A 



/3 2 A 2 ) 



/3A 
e{a 2 uj 2 



/3 2 A 2 ) 



(Bl) 



where := r)/[U 2 {n)(l — (n))] which vanishes in the 
half-filling particle-hole symmetric case. 

In this case, the bound states are still the solution to 
Eq. (|6j) if we replace a with a 

well as f3 with j3 

0(a 2 Lj 2 ~l3 2 A 2 ) 



2 (a 2 u, 2 



2 A 2) 



l-0 2 (a 2 w 2 -,3 2 A 2 ) i ' 

and also r\ with f] 



He 



~/3 2 A 2 ) ■ 



In the asymmetric case there are four solutions to 
the bound state energy Eq. ([6]). These are produced 



from a splitting in the two solutions obtained form the 
particle-hole symmetry. One can try a bound-state 
Ei, ~ \J 1 — t sin 2 cj) into the equation and come up with 
an effective transmission coefficient that depends on the 
Coulomb repulsion. 



/3 2 A 2 
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